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If) Using an SU(3) flavour symmetry breaking expansion in the quark mass, we determine the QCD 

*^ component of the neutron-proton, Sigma and Xi mass splittings of the baryon octet due to up- 

down (and strange) quark mass differences. Provided the average quark mass is kept constant, 
£SJ the expansion coefficients in our procedure can be determined from computationally cheaper 

simulations with mass degenerate sea quarks and partially quenched valence quarks. Full details 
and numerical results are given in [1]. 
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1. Introduction 

Isospin symmetry was introduced by Heisenberg in the 1930s to explain non-electromagnetic 
similarities between the proton and neutron. Nowadays, of course, this is ascribed to the u and 
d quarks having similar mass and the same strong - or QCD - interactions. This SU (2) flavour 
symmetry is not exact, there are isospin breaking effects, due to 

• the nid — m u quark mass difference which is a 'pure' QCD effect, 

• a QED component due to the different quark charges. 
As both effects are small then we can set 

M exp =M*+M QED , (1.1) 

where we denote by a * the 'pure' QCD component. There is an interplay between effects: electro- 
magnetic (EM) effects tend to make p heavier than n, but — m u works in the opposite direction 
and in fact dominates as the neutron is heavier than the proton, (M„ — M p ) exp = 1.293333(33) MeV, 
[2]. Including the s quark then the flavour symmetry group becomes SU (3) and the (pseudoscalar) 
mesons and baryons can be arranged in representations of this group. In Fig. 1 we show the lowest 
octet baryon and pseudoscalar states. States at the center, for example A(uds), L°(uds) have the 




Figure 1: The (lowest) octet baryon states, left panel and the octet pseudoscalar states, right panel. 

same quark content (quantum numbers) but different wavefunctions and can mix (if isospin is bro- 
ken) so we shall only consider states on the 'outer' ring here. As well as n - p mass splitting we 
now also have mass splittings involving the strange quark, (M^- — M J: +) exp = 8.079(76) MeV and 
(Ms- — M s o) exp = 6.85(21)MeV. These are all small differences (compared to the masses of the 
states) and the experimental precision is way beyond what we achieve here, but nevertheless we 
can qualitatively and reasonably quantitatively describe these splittings, as briefly described in the 
next section. For more details see [1]. 

2. Method 

The QCDSF-UKQCD strategy is to develop an SU (3) flavour symmetry breaking expansion, 
[3], from the flavour symmetric point down to the physical point. For the baryons on the outer ring 
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of the octet we have found up to NNLO 

M 2 (aab) = Mq +Ai(25jU a + 5jli fe ) +A 2 (SjU 6 - Sju a ) 
+ \B Q {8m 2 u + 8m 2 d + 8m]) 

+ B { (28 ii 2 + 8^) +B 2 {8pL 2 h - 8ii 2 )+B 3 (8ii b - 8^i a ) 2 
+ Co8m u 8m d 8 m s 

+ [Ci (25jU« + Sflb) + C 2 (5jU fe - 5jLl a )] (8m 2 u + + 

+ C 3 {8n a + 8pL b f + C A (8pL a + 8pL b ) 2 (8pL a - 8 ii b ) 

+ C 5 {8n a + 8pL b )(8pL a - 8n b ) 2 +C 6 {8n a - 8n b ) 3 , (2.1) 

while for the pseudoscalar meson octet, 

M 2 iab) = Ml n + a{8pL a + 8ix b ) 

+ A)K 5m « + 5m l + 8m D + Pi ( 8 ^a + +Pl(8lla ~ 8ll b ) 2 

+ Yo8m u 8m d 8m s + ji (8 \i a + 8\i b ){8m 2 u + + 8m 2 ) 

+ Yi{8pL a + 5jUfc) 3 + 73(5At fl + 8n b )(8n a - 8pL h ) 2 . (2.2) 

We have defined for the sea quarks 8m. q = m q — m with m = |(m M + m s ), where g G {a,b, . . .} 
so at the SU (3) flavour symmetric point 8m q = 0. From this definition this means that 8m s = 
—28m]. All the expansion coefficients are functions of m only. For the baryon or meson valence 
quarks we allow partially quenching, PQ, and set 8\l q = \L q — m (i. e. valence quark masses \L q ^ 
sea quark masses m q ). Of course on the unitary line when the sea and valence quark masses are the 
same then 8pL q — > 8m q . The quarks q = a, b, . . . are from {u,d,s}, so for example M(uud) = M p , 
M(dds) = M^- . We shall also need pseudoscalar mass results and the corresponding SU (3) flavour 
breaking expansion to determine the physical point: 8m* d , 8m* u and 8m* (where a * denotes the 
physical point). 

On the unitary line, singlet quantities have the property that the leading 0(8m q ) term vanishes. 
This allows a relatively simple definition of the scale, as practically we have shown, [3], that 
these quantities hardly vary in the interval from the flavour symmetric point down to the physical 
point. There are many possibilities for example for octet baryons (which are all stable under strong 
interactions), we may consider the 'centre of mass' of the octet 

X N = K m2 p +M 2 +M 2 + +M 2 - +M| +Ml-) =M 2 + 0(8m 2 q ) , (2.3) 

and similarly for the octet of pseudoscalar mesons 

X 2 = i(M| + +M 2 Ka +M 2 %+ +M 2 n - +Mi« +M 2 -) = M 2 K + 0(8m 2 q ) . (2.4) 

Using this we can form ratios M = M/X$ for S = N, % with expansion coefficients A,- =A,/Mq , a = 
cc/Mqx , ... for the SU (3) flavour breaking expansions. 

Note that as the coefficients of the SU (3) flavour breaking expansions are just functions of m 
alone, so provided m remains constant, the coefficients can be determined by ny = 2 + 1 simula- 
tions, when 8m u = 8m d = 8m\ rather than more expensive «/ = 1 + 1 + 1 simulations. Also an 
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additional advantage is that computationally cheaper PQ results can help to determine the coeffi- 
cients. 

Using 0(fl , )-improved clover fermions, [4], at /3 = 5.50 we have determined the appropriate 
point on the SU (3) flavour symmetric line (for the path to the physical point) and then used this 
point for PQ determinations of heavier baryon and meson masses (so that m obviously remains 
constant). Fitting these masses (and also including unitary data at the same constant m) then allows 
determinations of the expansion coefficients. In general these fits are functions of two quark masses 
8\l a and 8jJ.b- To avoid a 3-dimensional plot, as an illustration of these fits in Fig. 2 we show 



14 n . i ■ ■ ■ ■ i 70 
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Figure 2: The mass degenerate results, octet baryon M 2 (aaa) left panel and octet pseudoscalar meson 
M 2 (aa) right panel (using 32 3 x 64 sized lattices), both graphs versus the PQ quark mass 8fJ, a . Also shown 
are the functions given in eq. (2.5). 

(M 2 (aaa) - l)/(35jU a ) (left panel) and (M 2 (aa) - l)/(25jU a ) (right panel), together with the fit 
functions derived from eqs. (2.1) - (2.2) by taking completely degenerate quark masses 

M 2 (aaa)-\ ~ ~ „ „ ~ „ , M 2 (da)-\ ~ „ „ , 

=A 1 +B 1 5Ai a + fC 3 5/i a 2 , — \^ = & 1 +p 1 8Ha+4?i§ti- (2-5) 

As another example, we also consider the unitary results (i.e. 8pL q — > 8m q ) from the SU (3) flavour 
symmetric point down to the physical point. In Fig. 3 we show the baryon and pseudoscalar octet 
'fan' plots, where M N = M(lll), M z = M(lls), M z = M(lss), M Ns = M(sss) and M n = M{ll), 
Mk = M(ls), M rts = M(ss). The states Mjv f , M 7Jj are either not in the octet or are a PQ state and are 
not physical, but nevertheless can be used to help determine the expansion coefficients. The vertical 
lines are the rif = 2 + 1 pure QCD physical point, with the opaque circles being the determined 
pure QCD hadron mass ratios for 2 + 1 quark flavours. For comparison, the stars represent the 
average of the squared masses of the appropriate particle on the outer ring of the baryon octet, 
Fig. 1, i.e. M* N 2 {lll) = (MT p2 (ddu) +Mp Xp2 \uud)) /2, M* 2 (//s) = (M^ p2 (dds) +M^ 2 (uus))/2, 
M^ 2 (ssl) = (Mt xp2 (ssd) +Mt xp (ssu))/2. One immediate observation of Fig. 3 is that there is 
hardly any curvature in the data and that the NLO (i.e. quadratic terms) are sufficient. 

Finally note that the ^-scales used in Fig. 2 and Fig. 3 are very different \ 8mi\ ~ 0.01 <^L8jl a ~ 
0.5. (Indeed 8pL a ~ 0.5 is roughly at the charm quark mass.) The ability to use a large range for 
the PQ fits enables a much better determination of the fit coefficients (in particular the NLO terms, 
which are poorly determined in the narrow range \8mi\ ~ 0.01). 



4 



Isospin breaking 



R. Horsley 



1.8 
1.7 
1.6 
1.5 
1.4 
S 1-3 

o 

° 1-2 
/ 1.1 

O 

5 1.0 
0.9 
0.8 
0.7 
0.6 



-0.010 



-0.005 

8m, 



■ pure QCD 
A N(lll) 
-.Ellis) 
► S(SSI) 
TN.(sss) 
• sym. pt. 



o.ooo 



3.0 



2.5 



2.0 



fci 1.5 



0.5 



0.0 



-0.010 



-0.005 

5m, 



mil) 

<K(ls) 
Tt,,(ss) 
• sym. pt. 



o.ooo 



Figure 3: The baryon octet 'fan' plot, = Mjy /X^ (No = N, E, S, iV s ) left panel and the pseudoscalar 
meson octet 'fan' plot M^ g = M\ jX^ (%o = 7T> K, 7]. s ) right panel, both graphs versus The filled 

symbols represent mass values from 32 3 x 64 sized lattices while the opaque symbols are from 24 3 x 48 
sized lattices (and not used in the fits here). The common symmetric point is the filled circle. The stars (on 
the vertical line) are the estimated nf = 2 + 1 'pure' QCD physical points. The fits are eqs. (2.1) - (2.2). 



Of course, we are interested in mass differences here. So for the LO and NLO terms in eq. (2. 1) 
we have along the unitary line 

M n —M p = Middu) -M(uud) 

= (8m d - 8m u ) [A\ - 2A' 2 + (S\ - 2B' 2 )(8m d + 8m u )] , 
M £ - -M E + = M(dds) -M(uus) 

= (8m d - 8m u ) [2A\ -A' 2 + (2B\ -B' 2 + 3$' 3 )(8m d + 8m u )} , 
M s - — M z o = M(ssd) —M{ssu) 

= (8m d - 8m u ) [A\ +A' 2 + (S\ +B' 2 + 3B' 3 ){8m d + 8m u )] , (2.6) 

where the prime coefficients are simply related to the unprimed ones, [1]. Similarly we may invert 
the meson pseudoscalar expansion, eq. (2.2), to give 



8m d — 8m u 



K°~ M k+ 2(/3 1 +3fc,, 1 .~ 2 



a 



1+ 3«2 (^(%+%)-^ 



8m d + 8m u 



3a 



(l(M 2 K o+M 2 K+ )-Ml + ), 



(2.7) 



and then substitute in the baryon expansion eq. (2.6) to give the 'pure' QCD result. 



3. Results 

Performing this substitution gives the numerical results, [1], 

M n -M p = 0.0789(41)(34) (M 2 k0 -M 2 k+ ) [l +0.0817(92) Q(M|„ +m£ + ) -M* + )] , 
M L --M Z+ = 0.2243 (35) (92) (M 2 K0 -M 2 K+ ) [l +0.0077(30) (i(M|o +M 2 K +) -M|+)] , 
M E --M z o = 0.1455(24)(59) (M 2 k0 -M 2 k+ ) [l - 0.0324(50) (§(#£<, +M 2 k+ ) -M^+)] ,(3.1) 
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where M = M/X s , S = N, %. We see that the NLO corrections are small from +10% ~ —5% 
indicating that the SU (3) flavour symmetry breaking expansion appears to be a highly convergent 
series. (In eq. (3.1) the first error is statistical, the other is the total systematic error.) 

Note that eq. (3.1) is a 'pure' QCD result, and is the main result of this talk. We must now 
discuss what are the 'pure' QCD values of MjL, M% + , M 2 + . We know that EM effects are compa- 
rable to effects due to u - d quark mass differences. Dashen's theorem states that EM effects for 
charged mesons K + , n + are the same and for neutral mesons 7T°, K° vanish, Thus we can write 
Af^ 2 =M*J + fly, 2 = M*J « M* I , M e K x ? 2 = M* \ + ju y and M e ^ 2 = M* K 2 where a * denotes 
the 'pure' QCD 'physical' value, or M* K 2 -M* K l = {M 2 r0 -M| + ) exp + (1 + e y ) {M 2 n+ -M 2 ^)^ , 
where violations to Dashen's theorem (By = 0) are given by a non-zero Sy. We shall regard Sy here 
as a possible further systematic error, a typical value for it being Sy = 0.7, [5], giving about a 17% 
additional systematic error. 

Let us first investigate QED effects. From eq. (3.1), together with the numerical pseudoscalar 
meson masses from the previous paragraph we can determine the 'pure' QCD values. Then as we 
know the experimental values, [2] (as also given on page 2) then from eq. (1.1) we can determine 
the QED contribution to the mass splittings. This is shown in the left panel of Fig. 4. Thus this 
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Figure 4: Left panel: QED contribution to the baryon octet mass splittings, filled circles. Right panel: 
comparison of the n - p mass difference of the present result (QCDSF-UKQCD or bottom number) with 
NPLQCD, Blum et al., and RM123 [6, 7, 8] respectively (top to bottom). The filled circles use the QED 
determination of [9], while the filled square includes the full determination from Blum et al. The vertical 
dashed line is the experimental result. 

indicates that EM effects have the pattern 

n{ddu) — p(uud) ~ E (ssd) — E~(ssu) < 0, £~ (dds) — S + (uus) ~ . (3.2) 

Alternatively [9] gives a determination of electromagnetic effects of n - p of — 1.30(47) MeV 
(to be compared with —1.84(57) MeV here). In the right panel of Fig. 4 we compare our n - p 
mass difference including this determination of the QED contribution, (M„ — M p )* +QED , bottom 
result with the results of [6, 7, 8] (top to bottom). The filled square includes the full determination 
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from that reference. Despite the fact that QED effects are treated slightly differently in each work 
good agreement amongst the various determinations and with the experimental result is found. 

4. Conclusions 

We have introduced a method here to determine 'pure' QCD isospin effects in 

n-p, L--L+, S--S , (4.1) 

due to differences in u - d quark masses. This method involves developing a SU (3) flavour sym- 
metry breaking expansion keeping the average quark mass m constant. Advantages include the 
ability to use 2+1 simulations, i.e. m u =m c i = mi and use of computationally cheap PQ results. 
This expansion appears to be highly convergent, giving encouraging first results. Clearly the largest 
errors are due to unknown QED effects. For more details and numerical results see [1], 
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